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Abstract. - We study the electronic density of states in a mesoscopic superconductor near 
a transparent interface with a ferromagnetic metal. In our tunnel spectroscopy experiment, 
a substantial density of states is observed at sub-gap energies close to a ferromagnet. We 
compare our data with detailed calculations based on the Usadel equation, where the effect 
of the ferromagnet is treated as an effective boundary condition. We achieve an excellent 
agreement with theory when non-ideal quality of the interface is taken into account. 



Proximity effect (PE) is a phenomenon where the order parameter of a material in an 
ordered phase leaks to another material which has no such order. A number of experiments 
have investigated the PE in superconducting - normal (S-N) mesoscopic struct ures [p]|. Her e a 
finite pairing amplitude is formed in the N side, with the decay length = ^JhDph^k^ set 
by the dephasing of correlations between electrons and Andreev-reflected holes having slightly 
different momenta. At typical experimental conditions at low temperatures, the effect can 
survive up to large distances of several hundred nanometers. The experiments in non-magnetic 
metals have been generally well understood within the framework of the quasiclassical theory of 
inhomogeneous superconductivity 1^ . Another example of a proximity effect occurs in systems 
composed of ferromagnetic (F) and nonmagnetic materials, where the spin polarization leaks 
from F 1^ , but its decay length is typically very small. Going further, one may also consider the 
combination of two differently ordered materials, such as superconductors and ferromagnets. 

In this Letter, we study the mutual proximity effects of superconducting and ferromag- 
netic materials by measuring the local density of states in the superconducting side of the SF 
interface formed by superconducting Al and ferromagnetic Ni. Since the ferromagnetic and 
singlet superconducting order parameters try to exclude each other, the direct proximity effect 
is expected to get suppressed |Q] such that the pairing amplitude diffusing into a ferromagnet 
should decay over microscopic distances S,f = \/?iI?/27rfcB7curic due to the considerable dif- 
ference in the momenta of spin-up and spin-down quasiparticles. A spectroscopy experiment 
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in a low-Tcurio ferromagnet agreed with the theoretical picture of this short-range effect. 
However, some recent experiments |^-^ appear to indicate the existence of a long-range (i.e. 
a penetration depth comparable to ^^v) proximity effect into a ferromagnet and the situation 
has thus remained elusive. Here we report the first measurement on an inverse phenomenon, 
the modification of the BCS density of states in mesoscopic superconducting strips of Al under 
the influence of the proximity effect of a classical ferromagnet (Ni). This is rather a result 
from the expected suppression of the proximity effect: since the pairing amplitude vanishes 
at the SF interface, its value within a short distance from the interface is also much smaller 
than in the corresponding case of a finite proximity effect into a normal metal. 

We use tunnel spectroscopy with tunnel probes at fixed positions to measure the differential 
conductance on the S side of the interface. At zero temperature, and for a vanishing charging 
energy of the tunnel contact, the differential conductance G{E) with an N metal probe is 
proportional to the density of states (DOS). This way we can probe how the DOS in a 
superconductor is affected by the proximity effect. The DOS on the N side of an SN proximity 
system has been measured with fixed tunnel probes by Gueron et al. , who observed a dip 
in the N metal DOS up to one /im from the interface, with a quantitative agreement to theory. 
Several STM experiments capable of spatial mapping of the DOS have been carried out in 
short-coherence length layered superconductors |pl] -|l3| , and recently in Nb- Au structure . 

We used a sample structure shown in Fig. [l[ In the sample, we have two tunnel spec- 
troscopy probes: one at a distance of approximately 180 nm from the transparent SF interface 
(called the interface junction), and one 10 /im from the interface (called the bulk junction). 
The latter tunnel junction is assumed to be located at a site with properties of bulk Al. The 
sample was fabricated on oxidized Si wafers by electron beam lithography and shadow evapo- 
ration through PMMA/MMA copolymer mask. We first evaporated 20 nm of Al. To make the 
clean SF interface, we evaporated a 25 nm layer of Ni at a different angle. The depositions of 
the first two layers were done in a single UHV cycle at pressures below 5 • 10~^ mBar, avoiding 
any delay between the successive layers to make a transparent interface. The Al wires were 
then oxidized in 0.3 mBar of O2 for 5 minutes in order to make the opaque barriers for tunnel 
spectroscopy. Lastly, the tunnel probes were made by depositing 40 nm of Cu at a third angle. 
The process yielded resistances of approximately 1.5 MSI for both tunnel junctions. 

For electric transport measurements, the sample was cooled down to approximately 100 
mK in a plastic dilution refrigerator. Conductance measurements for the tunnel probes were 
performed through carefully filtered coaxial leads. To measure differential conductance, we 
used lock-in techniques with 20-30 /iV ac excitation added to a constant dc voltage bias applied 
to the tunnel probe to be measured. 

The experimental results are shown in Fig. ^. For the bulk junction, we obtained a BCS- 
like differential conductance, although the BCS peaks were somewhat lower than what would 
be expected simply by thermal broadening. The magnitude of the peaks was accurately 
accounted for only when we took into account the Coulomb charging effects (see below). For 
the interface junction, the peaks were lower, with a clear sub-gap conductance. 

In the normal state, superconductivity being suppressed by perpendicular magnetic field, 
we observed a conductance dip around small bias due to Coulomb charging effects. Both 
tunnel probes showed a dip of nearly 40 % due to a resistive environment created by the 
relatively resistive nickel wires. 

Diffusion constants for the wires or the interface resistance were not directly measurable in 
the sample. We used test samples, made in a similar process, to measure resistivities for thin 
wires of Al, Ni and Cu, with the results pAi = 1-7 fiilcm, pcu = 2.0 fiilcm and p^i — 32 fiftcm. 

The density of states near a SF interface can be calculated from the Usadel equation for 
the quasiclassical Green's functions in the diffusive limit. These functions can be parametrized 
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using the 0-parametrization j^. In the absence of supercurrents the retarded Green's function 
is & = cosh(0)f3 + i sinh(6')'r2, where are the Pauh matrices in Nambu space. In the case 
of translational invariance in the transverse directions, the Usadel equation then takes the 
form 

TiDdle = ~2iEsmh{e{E, x)) + 2iA(a;) cosh(6'(i;, a;)). (1) 

Here D is the diffusion constant and A(a;) is the superconducting pair potential, calculated 
self-consistently from 

A(a;) = A7Vo27i-A:BT^sinh6'(w„,a;), (2) 

where A^o is the density of states in the normal state, and the sum goes over the discrete 
Matsubara frequencies aj„ — {2n + l)nkBT and is cut off at the Debye energy. The attractive 
interaction between the electrons is characterized by the coupling parameter A, which is a 
nonzero constant for superconducting materials and zero otherwise. The parameter 9{ujn,x) 
is calculated from Eq. where the energy E is replaced by itUn- 

To account for a possible non-ideality of the interface between the two materials in the 
experiment, Eqs. (|l|,H) have to be complemented by a boundary condition recently derived by 
Nazarov relating the Green's function in the left side of the interface (Gi = G{x = 0^)) 
to the function in the right side (G2 — G{x — 0+)), 

p2 2%i G2,Gi 

aiGi9,G(0-) = a2G2a,G(0+) = ^V] ^ ^ ^. (3) 

„ 4 + r„(|G2,Gi}-2j 

In a junction with many transmission channels, the distribution of the transmission eigenvalues 
7n can be obtained from random matrix theory |l7| . In the case of a dirty interface it is given 

by @ 

where is the normal- state resistance of the interface per unit area. Integrating Eq. 
over the distribution of eigenvalues and over the cross sections of the wires, we finally get the 
desired boundary condition 

^ ' Rii ^ ' n v/i + cosh(0(O+)-0(O-)) ^ ' 

where = Rb/R(^ is the ratio between the normal-state resistance Rb = Pb^b of the 
interface with the cross section As, and the normal-state resistance R^ ~ ^s/ctnA of a piece 
of the right-hand-side (S) wire with length ^5, conductivity ajy and cross section A. The 
corresponding resistance of the wire in the left (F) side of the interface is denoted by Ri^ . 
From the solution 9{E, x) to Eqs. we obtain the local density of states through 

N{E, x) = NoRe{cosh{9{E, x))}. (6) 

In the ferromagnet, we assume the proximity effect to vanish within a few nm from the 
interface, and therefore set 9{x = 0^) = 0. In Fig. ^ we plot the DOS at varying values of 
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the interface parameter and at varying distances to the interface. Since couphng to the 
interface decreases with increasing rf,, DOS is very sensitive to the interface transparency. 

To calculate differential conductance of a tunnel probe, we have to take into account 
Coulomb charging effects due to the smallness of the tunnel junction capacitance. The charg- 
ing effects depend on the impedance of the electromagnetic environment in terms of the 
function P{E) which is the probability of exchange of energy between the tunneling electron 
and the environment. Forward tunneling rate through a single junction having the tunnel 
resistance Rt is given by (l9l 



nV)^4fr <^EdE'^^^^^^^§j^f{E)[l~ f{E' + eV)]P{E-E'), (7) 

where iVo is DOS of the tunnel probe, taken as constant, and Ns is the proximity-affected 
DOS on the S side. The function P{E) is calculated from an integral equation assuming a 
purely resistive environment ]2C[ |. 

For both tunnel probes, the parameters a = Rg/Rcnv and the charging energy of tunnel 
junction capacitor Ec were determined as the values that gave the best fit to the normal- 
state conductance. Here, Rq = /i/4e^ ~ 6.5 hfl is the superconducting resistance quantum, 
and i?env is the resistance of the electromagnetic environment seen by the tunnel junction. 
Normal-state data of both junctions were fitted with the same values. These values were then 
used to fit A in bulk Al data, and we obtained a faultless fit with A ~ 0.22 mV and with the 
temperature T = 100 mK that was recorded in the experiment. 

The distance between the interface and the first probe had some ambiguity of ±50 nm 
because the tunnel probes had a finite width of 80 nm, and because the two metal films at 
the interface overlapped each other by approximately 50 nm. Average of the DOS over these 
distances was practically equal to the DOS in the middle of the tunnel probe, and we used 
the average distance (this distance is about 1.2 ^s) in the analysis. 

The normal-state and bulk-junction measurements yielded all the other important param- 
eters of the system except the interface parameter rb. For example, the electronic temperature 
in both bulk and interface junctions can be assumed the same due to the similarity of the probe 
resistances. Typical resistances reported for SF interfaces similar to our sample have been in 
the range 2017 — 3017 For our sample made in UHV conditions, we expect Rb — 10 fl. 
Since i?^ « 2 il, we would have r?, ~ 5. 

With a reasonable value of the interface parameter r;, = 3, we get an excellent fit with the 
theory, as seen in Fig. I Note that this fit was obtained without taking into account such 
effects as the magnetic field due to the ferromagnet, penetration of the spin polarization in 
the superconducting side, or the emergence of a triplet proximity effect pl| . 

The effect of a stray magnetic field from the ferromagnet can be examined as follows. We 
expect the size of magnetic domains in the film to equal the thickness of the film, d = 25 
nm. Since the dipolar field from a single domain decays with distance r as (d/r)^, at relevant 
distances r ~ 200 nm, the saturation field of Ni of 0.6 T has dropped down to 1 mT. Compared 
to the much larger critical field of the films (70 mT), this is unlikely to explain the enhancement 
of DOS at low energies. To gain further confidence that the observed DOS was not affected 
by the stray magnetic field, we measured the conductance at various externally applied field 
strengths. At fields up to ~ 20 mT, the zero-bias conductance stayed practically constant 
while the width of the gap decreased. This behavior was similar in both the interface junction 
and the junction at bulk Al. Thus, stray fields even significantly larger than the estimated 
would not be likely to affect the zero-bias DOS, but would only reduce the gap. The case is 
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strengthened by the fact that the gap was not reduced more than expected from the Usadel 
theory. These arguments also favor the idea that the domains were indeed of the estimated 
size. 

In conclusion, we have measured the inverse proximity effect in a mesoscopic superconduc- 
tor in contact with Ni, and compared the data with a self-consistent calculation based on the 
Usadcl equation where we assume Ni to suppress the PE in a microscopic scale. Wc obtain a 
good agreement for the electronic density of states close to the interface only by taking into 
account the effect of the non-ideal interface. 
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Fig. 2 - Normalized differential conductance G/Go of the SIN junction, when properties of S are 

affected by the proximity of a ferromagnet. The circles are the experimental data at a distance of 
180 nm (1.2 ^s) from the interface, and that for the bulk Al, shown for reference (measured at T = 
100 mK). The data are normalized to the conductance value Go measured at 1 mV. The solid curves 
are the best fits to theory, with a = 7.5 (implying iienv = 870 O), Ec/A = 1.5, ri = 3. The gap of 
Al was A = 0.22 mV. 
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Fig. 3 - DOS of SF structures calculated on the S side, (a) with an ideal interface, ri, = 0, at 
distances 0.5 ^§...3 ^5 from the interface; (b) at a fixed distance x = ^5 with varying interface 

parameters r(, = 1...3. In the calculations, the S order parameter was assumed not to penetrate F, 
i.e. DOS in F was taken as constant A^o. The inset shows the self-consistent pair potential. 



